Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 2 February 2008 (MN KTeX style file v2.2) 



Particle Energies and Filling Fractions of Radio Bubbles in Cluster 
Cores 

R.J.H. Dunn * and A.C. Fabian 

Institute of Astronomy, Madingley Road, Cambridge CB3 OHA 



o 
o 

(N 

CD 



> 

in 
vn 
o 
o\ 
o 

o 



Oh 
I 

o 



X 



2 February 2008 



ABSTRACT 

Using Chandra images of cluster cores with clear radio bubbles, we have determined k, which 
is the ratio of the total particle energy to that of the electrons radiating between 10 MHz and 
10 GHz. Radiative and dynamical constraints on the bubbles indicate that the ratio of the 
energy factor, k, to the volume filling factor, /, lies within the range I k/ f < 1000. 
Assuming pressure equilibrium between the radio-emitting plasma and the surrounding X- 
ray gas, none of the lobes have equipartition between relativistic particles and magnetic field. 
There is no evidence for any dependence of the upper limit of the k/ f ratio on any physical 
parameter of the cluster or the radio source. The distribution of the upper limit on fc// appears 
to be bimodal, the value for some clusters being ^ 3 and for the others ^ 300. We show that 
this is may due to the composition of the jet which forms the bubbles, the variation in the 
volume filUng fraction or variation in the amount of re-acceleration occurring in the bubble. 
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1 INTRODUCTION 

Radio lobes which emit synchrotron radiation contain relativistic 
electrons and magnetic fields. Disentangling the pressures of each 
component is difficult, and traditionally it has been ass umed that 
there is equipartition between the particles and the field iBurbidgj 
Il959h . which corresponds closely to the minimum energy condi- 
tion. Now, however, when the bubbles are embedded in an Intra 
Cluster Medium (ICM), this degeneracy can be removed by mea- 
suring the thermal pressure of the ICM and assuming pressure 
equilibrium between the bubbles and the surrounding X-ray gas 
as well as the lack of strong shocks. In particular, FR I sources 
in low redshift clusters imaged with Chandra often show holes 
in t he X-ray emissio n coinciding wit h the radio lobes (e .g. Hydra 
A. iMcNamara et al. 2000); Perseus, jpabian et alJl2000l) : A2052, 
j Blantonet jdlE oOD: A2199, (Johnstone et alj|2002h : Centaurus. 
iSan ders & Fabian 2002)), the first of which was discovered in the 
Perseus cluster with ROSAT ( Bo hringe r et al. 1993). A recent com- 
pilation is given bvlBu-zanetalJ l2004) . 

Here we perform a detailed study of a sample of clusters which 
have clear radio bu bbles. Following the approach first detailed in 
iFabian et al. ( 2002). we determine fc//, where fc is the ratio of the 
total relativistic particle energy to that in electrons emitting syn- 
chrotron radiation between 10 MHz and 10 GHz, and / is the vol- 
ume filling factor of the relativistic plasma in the bubble. 

Fabian et al. (2002) calculated fc// for the northern radio bub- 
ble in the Perseus cluster and obtained a value of 180 < fc// < 
500. We find a value that is in broad agreement with theirs, and find 
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similar values for M84 and PKS 1404-267. However the remaining 
clusters in our sample, including A2052, A2199, A4059, Centau- 
rus and Hydra A, appear to have much lower values, fc// appears 
to be either in the range ~ 1 — 10 or 100 — 1000. We discuss 
the implications of this result on the evolution and formation of the 
bubbles. We use Hq = 70 km Mpc^^ throughout. 



2 DATA ANALYSIS 

We use standard synchrotron theory to quantify the properties of the 
particles present in the bubbles, following the analysis presented 
in Fabian et al. (2002). We give below some of the formulae re- 
quired in order to define the variables. The total energy in relativis- 
tic electrons radiating between u\ to V2, with a spectral index a 
{S{u) oc v"), in a magnetic field B, producing a flux density Sv at 
u, is 



An X 10' 



12/ cz 



1 + 



2 



,0.5+a 



-B 



-3/2 



.HoJ \ 2/1/" a + 0.5 
^ aB~^'^ erg, (1) 

where Ho = 70kms~^ Mpc~^. It has been assumed that the 
particle energy distribution extends from = 10 MHz to U2 = 
10 GHz as the radio images analysed were taken at frequencies 
between 1 and 8 GHz. Therefore the total energy in particles and 
magnetic field is 

o2 

Etot = fc£c + Vf— = afcB-=*/' + 6/B' erg, 
Stt 

where V is the volume of the bubble, and / represents the volume 
filling factor of the relativistic plasma, k accounts for the additional 
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energy from relativistic particles accompanying the electrons that 
radiate above 10 MHz and any non-relativistic component (k = 1 
for an electron-positron plasma emitting only in the above wave- 
band; a typical value used in the literature is k — 100). The bub- 
bles are allowed to be aspherical, where ri ("length") is the radius 
along the jet direction and r„ ("width") is the radius across the jet 
direction, where a symmetry axis has been assumed along the jet, 
hence the volume is \/ = 47rrir^/3. 

With the condition that there is simple equipartition between 
the energy present in particles and that present in the magnetic field, 
the magnetic field strength is 

For the minimum energy magnetic field, replace with 
(3a/46)^'^^ in the above equation. 

If the relativistic gas is in equilibrium at equipartition with the 
thermal pressure from the gas in the rims surrounding the radio 
lobes (Pth), then the equipartition value, fc//cq, can be obtained. 

From now on we assume that the plasma is not in equipartition 
but still in pressure equilibrium, hence k/ f can be found from 

/ V 8tt / a 

To calculate k/f from Equation |2| estimates of the magnetic field 
strength in the bubbles are obtained from radiative and dynamical 
constraints. 

A minimum value of fc// — 1 can be understood for an 
electron-positron plasma which fills all of the bubble. However, 
an extremal value can be determined by differentiating Equation 
|2| and finding the corresponding magnetic field gives a maximum 
value, fc//max, which is 50% greater than the equipartition value. 
The absolute maximum value for the magnetic field is when k/f — 
1 from when (Pth ~ -B^/Stt) ~ 0, which is 1.53 times greater than 
the magnetic field at fc/Jmax, and 1.15 times the equipartition mag- 
netic field. The field for k / /max is the limit up to which the term 
can be ignored. Any further increase in B and this term becomes 
dominant and k/f decreases until it equals zero. After that the mag- 
netic pressure is such that the bubble would be over pressured, even 
with k — 1, and for pressure equilibrium fewer particles than ob- 
served would be required. 

A limit on the magnetic field present in the bubbles can be de- 
duced from the fact that GHz radio emission is seen throughout the 
bubbles. Therefore the synchrotron cooling time of the relativistic 
electrons present is 

Uync = 2.7 X lO^BZTu^^^^ yr, (3) 

where B — 10~^B_5 G and ly = 10^ ug Hz, has to be greater than 
the age of the bubble, assuming that there is no re-acceleration of 
the electrons in the bubble (for further discussion see Section lOt . 
However Equation|3|is only valid for strong enough magnetic fields 
(for further discussion see Section lOt . 

The age of the bubble can be determined from the fact that the 
X-ray rims surrounding the radio bubbles appear not to be strongly 
shocked, hence they must be expanding at less than the sound speed 
of the gas in the rims, Cs = \/ jkT/ fimn, where 7 is the ratio 
of the heat capacities and fi = 0.62. And therefore the age of the 
bubbles must be greater than fsound ~ 2ri/cs where ri is the radius 
of the bubble. Twice the radius has been used as the bubble is blown 
from the radio source at the centre of the galaxy, which is at one 
edge of the bubble, rather than the centre, and so the "front edge" 



travels twice the bubbles' radius. This edge has travelled furthest 
during the lifetime of the bubble and so provides a limit on its age. 

If the bubbles rise upwards at their buoyancy velocity, = 
^2gV/SCu, where S is the cross-sectional area of the bubble, 
V is the volume, g — GM{< -Rdist)/iidist for the bubble (cen- 
tre) being at Rd ist from the cluster co re and Cd = 0.75 is the 
drag coefficient JChurazov etai]|200ll) . then the age of the bub- 
ble can be estimated as ibuoy = -Rdist/i'b, the travel time to their 
current position. The enclosed mass was estimated from cluster 
mass profiles or using a linear interpolati on from the Abell ra- 
dius and mass as listed in Reiorich & Bo hringeJ i2002h . The lat- 
ter is probably only correct to a factor of a few, though the un- 
certainty was taken into account in the calculation. The age of the 
bubble can also be estimated from the time requi red to refill the dis - 
placed volume as the bubble rises upward tMcNamara et alj200(t) . 
t.cflii = 2Rii,ty^r/GM{< i?dist). 

As the synchrotron lifetime of the GHz electrons must be 
greater than the age of the bubble, limits can be placed on the mag- 
netic field, assuming that ug — 1. 

The value obtained shows whether the equipartition solution 
is possible and allows k/f to be calculated from Equation |2| The 
limits obtained for k/f may be higher than the one predicted for 
equipartition. However, the equipartition con dition is determined 
from the magnetic field strength, and Fig 7. of iFabian et alj J2b02h 
shows that it is possible that limits on k/f are larger than the 
equipartition value even though the magnetic field is less than its 
equipartition value. 

In the Perseus cluster there are weak shocks visible in the X- 
ray images 1 Fabian et al. 2000). Assuming that these have been cre- 
ated by the earlier supersonic expansion of the bubble, then a max- 
imum age for the bubble can be estimated. Currently these shocks 
are travelling at the local sound speed, however, in the past the bub- 
ble that caused them was travelling supersonically, and so the av- 
erage speed of the shocks over the bubble lifetime will be greater 
than the sound speed. Hence the age of the bubble must be less than 
^shock = -Rshock/cs where i?shock is the shock radius, centred on 
the central radio source. 

Although this is an upper limit on the age of the bubble, there 
is no requirement that the synchrotron lifetime of the electrons is 
less than this value. However, in the assumption that it is, then this 
upper age limit, gives a lower limit on the magnetic field, and hence 
a lower limit on the value for k/f. In the Perseus cluster, which is 
the only cluster in the sample with a shock front visible, this limit 
together with the value obtained from the condition that there are 
no strong shocks at the rims gives the lowest range of fc// that can 
be u nderstood from physi cal arguments. 

iFabian et alJ i2003lJl showed that some of the Hq filaments 
seen in the Perseus cluster could be interpreted as streamlines be- 
hind the western "ghost" bubble. This indicated that the flow is 
laminar, and implies a Reynolds number of less than 1000. Using 
this, a value of 4 x lO'^'^ cm"^ s^^ was obtained for the kinematic 
viscosity. The Reynolds number. Re, was calculated for each bub- 
ble analysed, assuming the viscosity was the above value, and also 
a lower bound on the viscosity was calculated from the limit on the 
Reynolds number of 1000 assuming that the flow is laminar in each 
cluster. 

2.1 Non-uniform Magnetic Fields 

It is unlikely that the magnetic fields present in the bubbles are com- 
pletely uniform over the entire bubble. The effect of non-uniform 
magnetic fields was investigated very simply, using a centrally 
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Table 1. : Effect of Non-uniform Magnetic Field on k/f 



Cluster Lobe 


Timescale 


Uniform Field 


Non-uniform Field 


Perseus North 


Sound 


573 


1083 




Shock 


390 


735 


Perseus South 


Sound 


1025 


1937 




Shock 


762 


1436 



Notes: 

The limits on k/f for the sound speed and shock calculations from the 
Perseus Cluster with a uniform and non-uniform magnetic field. 



peaked linear distribution (B oc — [rp. This gave limits on k/f 
that were a factor of two greater than if a uniform magnetic field 
was used, for the case of Perseus (Table0. It was checked that the 
calculation used for the non-uniform field produced the same lim- 
its on fe// for the a uniform field as were produced from the other 
calculations. 



2.2 Energetics 

The energies and powers were also calculated, using E = PV for 
the energy, and V — E/t for the powers using all of the timescales 
calculated above. However, if the energy required to expand the 
bubble as well as the energy in the bubble is accounted for, then, 
for slow expansion rates: 



E - 



7 - 



-PV + PdV 



7 



■PV, 



7 - 



where V is the volume of the bubble and 7 is the mean adiabatic 
index of the fluid in the bubble (5/3 for non-relativistic gas or 4/3 
for relativistic gas). For the relativistic case (assumed as there is 
synchrotron emission in the cavity) this increases the energy by 
a factor of 4 (5/2 for non-relativistic), and hence also the power. 
Therefore the values for the energy in Table|3|should be multiplied 
by the appropriate 7/(7 — 1). 



3 DYNAMICAL CONSTRAINTS 

PdV work is done by the bubble on its surroundings during its ex- 
pansion, equating to a power P ~ PthfV/t. The jet must be pow- 
erful enough to inflate the holes, and an estimation of this power 
can be obtained by considering the ev olution of the expanding bub - 
ble in a medium of constant pressure JChurazov et all200lLl200ol) . 
We modify their expressions for the expansion of the bubble to in- 
clude the effect of gas clumping through the volume filling factor. 
From the internal energy of the hole and the compression of the 
surrounding gas during subsonic expansion 



Vt 

Pth/ 



= V- 



7 



7 - 



(4) 



where 7 for the radio-emitting plasma is assumed to be 4/3. Sec- 
ondly, from the condition that there are no strong shocks at the 
bubble rims, the expansion rate must be subsonic, and so 



r 



< 367rc, 



7 



7 - 



= G. 



(5) 



Finally, assuming that the bubble has not detached from the source, 
which seems to be the case for the majority of the clusters, then its 




2x10 



'max IXlO' 

Time, (yr) 



2x10' 



Figure 1. The dynamical constraints given in Sectionl3lshown for the north- 
em Perseus bubble. The vertical lines show the minimum tmin, obtained 
from the combination of the energy and sound speed equations |4} and js}, 
and maximum tmax ages of the bubble, obtained from the combination of 
the energy and buoyancy equations and The inclined lines show the 
limits on the power of the source, the upper shallow line is from Equation 
l^and the steep line is from EquationQ The dotted shallow line represents 
the minimum power required to move the gas to make the hole (Pn^V/t). 
In all the equations, / has been taken as 1, and Pth is the value obtained 
from the temperature and density of the bubbles. 



velocity is less than the buoyancy velocity, giving 



(6) 



Combining these expressions produces other limits on the ages of 
the radio bubbles — J4j and Js} produce a minimum age for the 
hole, tniin = ^f"^ and ^ lead to a relation 



P < H 



-2/3p5/3 



Ptht- 



(7) 



There is also a minimum power for the source, from PthV/t for 
f — 1. The intersection of this with gives the maximum age of 
the hole as 



imax — 



P5/3H-2/3 



V 



These relations are only valid for subsonic expansion. For Perseus 
at least, there exists a weak shoc k indicating that the re has been su- 
personic expansion in the past i Fabian et all2003'3) . Therefore the 
agreement between the timescales given by these relations, with the 
ones drawn from the other physical arguments may give a crude in- 
dication of what proportion of its lifetime the bubble has spent in 
subsonic expansion. The allowed region for the bubble from these 
constraints in the Power-Age plot in Fig.Qis the highlighted poly- 
gon defined by the minimum and maximum ages, the power limits 
from Equations|4]andQand the minimum power for the source. 



4 RESULTS 

The individual source properties are listed in Table^l In A2199, the 
inner jet-like feature visible on the radio overlays in lohnstone et all 
(1998) has small X-ray depressions associated with the knots at the 
ends of the jet (Inner East (I-E) and Inner West (I-W)), and these 
were analysed separately from the extended primary emission. If 
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Figure 2. The value oik/ f calculated from the sound speed limit for each 
cluster analysed, along with the uncertainties arising from the uncertainties 
in a (the dotted bars) and from the uncertainties in the other physical pa- 
rameters of the source (the solid bai's). The dotted line shows the minimum 
value of fc// possible from the assumptions used in the calculations. 



Figure 3. The value of fc// for each cluster analysed. The I, A, • and (} 

symbols denote the k/ f values from the sound speed, shock, buoyancy and 
refilling timescales respectively. The X symbol denotes the equipartition 
values 



the source has prominent X-ray rims around the radio bubbles, e.g. 
Perseus, then these have been used to define the size of the bubbles. 
Otherwise the spatial extent of the radio emission has been used to 
define the size of the radio bubble. In some cases the GHz radio 
emission does not entirely fill the X-ray defined holes - e.g. A2052. 
In Cygnus A, the fluxes used in the calculation do not include the 
hotspots at the end of the jet, to reduce the effect of re-acceleration 
of electrons (for further discussion see Section lSHI . 

The limits on fc// calculated from the different timescales are 
tabulated in Table|3| Ignoring the uncertainties in the limits on fc//, 
Perseus has the smallest range in fc// as the two constraints come 
from the weak shock (fc//shock) and the lack of strong shocks in 
the rims (fc//sound)- However the ranges in fc// for the Northern 
(390 < fc// < 570) and Southern (760 < fc// < 1025) bub- 
bles do not overlap, which is surprising. However, if the uncertain- 
ties in the individual limits are also included, then the ranges are 
much larger, and do overlap (200 < fc// < 1020 and 390 < 
fc// < 1830). The values obtained are in agreement with t he range 
of 180 < fc// < 500 obtained by iFabian etall <2002t) for the 
northern bubble. This is assuming that the synchrotron lifetime is 
less than the maximum age of the bubble, but there is no physical 
requirement that this is the case. All other clusters have an upper 
limit on fc// arising from the sound speed limit and two other limits 
from physical arguments (buoyancy and refilling timescales). 

The maximum and minimum values in the table above come 
from the uncertainties in the values of the spectral index, a, over the 
radio lobe. The resultant uncertainty in fc// from the uncertainty in 
Q can be greater than an order of magnitude as the spectral index 
appears as an exponent in Equation Q There are uncertainties in 
fc// which come from the uncertainties in the other observed val- 
ues, however for most of the clusters these are much smaller than 
the uncertainties from the those in a. The limits on fc// calculated 
from the lack of strong shocks in the rims (sound speed limit) is 
plotted in Fig.|2|with the uncertainties in the limits on fc// aris- 
ing from the ranges in a and those from other uncertainties plotted 
separately. The dotted line demarcates the minimum value possible 
for the assumptions used in the calculations. Each lobe is plotted 
separately, identified by the same label as in Tables|2|and|3| 

The uncertainties in the spectral index are generous — a 



varies across the radio lobes, and so an average value has been 
taken, with an uncertainty to cover the range. In some cases the 
values stated by different authors do not match, and again an aver- 
age has been taken, with the range encompassing the values. 

The variation in the limits on fc// obtained from the differ- 
ent timescales is shown in Fig.|3| On the whole, the fc// limiting 
values from different timescales fall within the range of fc//sound 
given by the uncertainties in a. It was not expected that the val- 
ues would be exactly the same, but for about half of the clusters 
the differences are less than ~ x2.0. The remainders' are less than 
~ x8 and these disagreements arise primarily from differences in 
the timescales calculated (for further discussion see Section l4m . 
As the values obtained for the limits on fc// from these different 
arguments are on the whole similar, it implies that the actual limits 
on fc// are around those values. However all the values have uncer- 
tainties of the same size as for the sound speed (they have not been 
plotted for clarity), which reduces the reliability of this conclusion. 

For most of the bubbles, radio images at different frequencies 
were obtained, and a separate limit on fc// was calculated for each. 
The resultant values usually were all but identical, and certainly 
agreed to within the uncertainties arising from all values apart from 
the spectral index, however they have not been tabulated in this 
work for brevity and clarity. It was checked to see if the frequency 
at which the flux was measured had some effect on the limit on 
fc//, but none was found. In the clusters where two bubbles are 
visible, their fc// limits broadly agree to within the uncertainties 



4.1 Timescales 

The timescales (ages) for each lobe were calculated during the 
course of this work. The differences between our values and those 
presented in the literature come mainly from the different sizes of 
the bubbles used, and different cosmologies. They have been tabu- 
lated in the Appendix along with the powers of the sources which 
have been calculated from these timescales. 

For some of the bubbles the timescales obtained by the dif- 
ferent estimates differ by up to a factor of 10. These discrepancies 
are probably due to the assumptions used in the calculation of the 
timescales. 

The sound speed timescale relies on the fact that the bubble is 
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Table 2. : Source Properties 



Cluster 


Lobe(i) 


Redshift 


a 


log(P5GHz) 

(W/Hz) 


d(2) 
-^dist 
(kpc) 


n 

(kpc) 


(kpc) 


^shock 

(kpc) 


(IQI^Mq) 


kT 
(eV) 


(cm~^) 


References 


A133 


Rel, R 


0.054 


-2.0 ± 1.0 


23.0 


40.0 


50.0 


20.0 


- 


8.0 


2400 


0.01 


1,2,3,4,5,6,7 


A2052 


N, X 


0.035 


-1.6 ±0.3 


24.7 


10.0 


10.0 


10.0 


- 


4.0 


1100 


0.049 


7,8,9,10,11,12,13,14 




S, X 








15.0 


15.0 


15.0 


- 


5.0 


1100 


0.049 




A2199 


E,R 


0.031 


-1.9 ±0.3 


24.3 


20.2 


16.9 


12.2 


- 


1.6 


3100 


0.017 


13,15,16,17 




W,R 








21.5 


20.9 


14.4 


- 


1.7 


3200 


0.017 






I-E, R 








4.23 


1.0 


1.0 




0.08 


2050 


0.056 






I-W,R 








4.23 


1.5 


1.5 


- 


0.08 


2050 


0.056 




A4059 


N, R 


0.050 


— 1.4 ± 0.3 


24.0 


5.44 


5.44 


2.98 




1.5 


2000 


0.028 


7,13,14,18,19 




S,R 








3.92 


3.92 


4.48 


- 


1.1 


2000 


0.028 




V^^ilLuLll Lid 


E R 


0.010 


— 1.3 ± 0.2 


25.8 


1.45 


1.45 


1.13 




0.020 


800 


0.09 


13 20 




E,R 








2.00 


1.88 


1.02 


- 


0.025 


1000 


0.08 




Cygnus A 


E,R 


0.056 


-1.0 ±0.3 


27.7 


44.1 


30.3 


23.7 


_ 


3.0 


4000 


0.05 


13,21,22,23,24 




W,R 








48.0 


33.9 


25.2 




5.0 


5500 


0.04 




Hydra A 


N, R 


0.052 


-1.2 ±0.4 


26.3 


17.8 


20.0 


8.5 


- 


2.0 


2900 


0.041 


13,25,26,27 




S, R 








22.2 


16.7 


6.9 




2.0 


3000 


0.036 




M84 


N, R 


0.0035 


-0.51 ±0.02 


22.0 


2.7 


2.1 


2.5 


- 


0.2 


650 


0.04 


28,29,30 




S, R 








3.8 


4.0 


2.6 




0.2 


650 


0.04 




M87 


Bud, R 


0.0044 


-1.2±0.4 


24.4 


3.11 


1.31 


1.31 




0.14 


1500 


0.20 


13,31 


MKW3s 


Rel, R 


0.045 


-2.7 ±0.05 


23.8 


56.4 


12.9 


17.9 




16.0 


3200 


0.004 


7,32,33 


Perseus 


N, X 


0.018 


-1.0 ±0.2 


25.9 


8.15 


8.15 


8.15 


24.0 


0.24 


2604 


0.064 


13,34,35,36,37,38 




S, X 








8.89 


8.89 


8.89 


24.0 


0.26 


3175 


0.055 




PKS 1404 


E,X 


0.022 


-0.39 ±0.03 


23.8 


3.94 


2.29 


1.58 




0.40 


930 


0.056 


39,40,41,42 




w,x 








2.75 


1.33 


1.7 




0.40 


820 


0.066 
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iTerasranta et aL t2001.): (30'I.Cao & Rawlings i2004.) : r3n.Forman et al..t2003.): (32).Mazzotta et al..(.200Z): (■3 3 VMcNamara et al. (1990): (34) Fabian et ail 
l2002ft: (35) iFabian et alj l2OO0ft: (36) iFabian et alj l2003th: (37) IVoigt & Fabi'^ <2004 : (38) IPedlar et all <199cl) : (39) brinkwater et ai] <1997ft : (40) 
Ijohnstone et alj|l998ft : (41) lBrown & BumJ<199lh :" 42) lBettoni et alj|2003ft 
Notes: 

(1) The codes for the Lobes are N — Northern, S — Southern, E — Eastern, W — WesteiTi, X — sizes from X-ray image, R — sizes from Radio image, I-EAV — Inner 
lobes in A2199, Bud — as described in Forman et al. ( 2003), Rel — Relic source (for A133 as described in Fuiita et al. ( 2002)). 

(2) All the yalues giyen in the aboye table except the radio power haye an uncertainty associated with them. Except for the spectral index, they are not quoted 
as they haye limited effect on the calculated yalues. The effect of the uncertainties in a on the resultant uncertainties in A:// is large and so are stated here. For 
further discussion see text. 



not currently expanding supersonically and gives the time for the 
expansion of the bubble. From the observation that there are no 
strong shocks at the rims we have assumed that the average expan- 
sion speed was less than the current sound speed over the lifetime of 
the bubble. However, as we assume that the bubble is produced by a 
relativistic jet and the presence of the shock in the X-ray images of 
the Perseus cluster, the earlier expansion of the bubble must have 
been supersonic, and so this lower limit on the timescale may be 
an overestimate. The lac k of strong shocks in the rims still allows 
some supersonic motion. lsianton et alJi200l'l) state that a < 1.2 



for A2052, so the sound speed may be an underestimate of the ex- 
pansion speed. 

The buoyancy timescale used is the time taken for a bubble 
of the size observed to rise at its buoyancy velocity from the cen- 
tre of the gravitational potential to its current position, i.e. by its 
own radius. This is calculated assuming that the surrounding ICM 
is at a constant pressure and density. However as the bubble is ex- 
panding as it ages, the initial stages are probably supersonic, the 
medium through which it rises is likely to be non-uniform and in 
some cases the masses used in the calculation are estimates; the 
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Table 3. Physical fc// Values 



Cluster 


Lobe(i) 


Pressure 
(eV/cm^) 


Re 


Viscosity 
(10^^ cm^ s-i) 


Energy '3) 
(10^8 erg) 


k/ /oq 


J sound 


^/ /"shock 


^//buoyancy 


fe/ /rcflll 


A133 


Rel, R 


52.8 


4106 


16.4 


20.8 


90.6 


8.79 t,AO 




c:7 o 2760 
'^'■■^ 0.68 


17.5 8« 


A2052 


N, X 

s,x 


119 
119 


367 
551 


3.83 
4.69 


2.3 
7.9 


26.6 
81.3 


4.73^537 
9.68 




36 3 l^^ 
■^""^ 6.19 

79.7 1^6 


12.6 68J 

94 131 
4.10 


A2199 


E,R 
W,R 
I-E,R 
I-W,R 


116 
120 

253 
253 


754 
904 
50.1 
75.2 


3.77 
4.90 
0.08 
0.15 


4.2 
7.1 

0.005 
0.017 


21.5 
39.8 
2.11 
4.32 


n oc 14.0 

rr 74 21.0 
3. 1.50 
2.40 l^il 
CO 18.0 
0.68 


- 


7.90 28.8 
14.8 13.9 

0. 56 git 

1 40 7.06 

1. ^*^ 0.274 


3 02 11.0 
•^■^^ 0.79 

4 07 18.1 

1.30 

0.63 

1.06 5.32 


A4059 


N, R 
S,R 


123 
123 


148 
222 


0.69 
1.22 


0.065 
0.22 


9.97 
25.0 


A 1 18.0 
^■^^ 0.86 
14 1 61.1 
1^.1 2.94 


- 


14.9 i^i^ 
26.6 1% 


Qir 38.3 

91 c; 93.0 
4.47 


Centaurus 


E,R 
E,R 


158 
176 


35.3 
35.7 


0.16 
0.13 


0.0045 
0.0037 


1.57 
2.00 


1 91 3.16 
1.^1 0.44 
1 26 3.28 


- 


9 IS 5.67 
^.1° 0.79 

2.03 5.29 


0.86 2.25 

0.88 2.28 


Cygnus A 


E,R 
W,R 


440 
484 


1552 
2066 


6.21 
8.26 


116.0 
152.0 


35.5 
46.9 


1 50 486 

1 04 6.27 
1.^* 0.48 


- 


1 09 6.23 
1.^'^ 48 
9 70 8.84 
'^■'■^ 0.68 


0.84 2.72 

1 10 3.86 
l.l'^ 0.30 


Hydra A 


N, R 
S,R 


264 

238 


507 
416 


3.64 
2.16 


3.2 
1.5 


6.84 
4.45 


n 3.11 
'^■'^'^ 0.07 
n 47 2.65 
I'.*' 0.06 


- 


9 1 R 12.2 
^.1° 0.29 
n Q1 5.13 


0.80 4.53 

50 0.84 
'J.^'J 0.07 


M84 


N, R 
S,R 


57.2 
57.2 


71.8 
73.5 


0.52 
0.52 


0.019 
0.020 


205 
196 


201 207 

110 


- 


101 1'A 

298 i?i 


268 

202 208 


M87 


Bud. R 


660 


56.2 


0.22 


0.029 


56.2 


23.9 1^, 




19.7 lo^^ 


12.6 44.7 


MKW3s 


S,R 


28.2 


1118 


5.50 


3.2 


2.04 


1 00 1.66 
^"^^ 1.05 




0.78 g:i 


0.76 g:if 


Perseus 


N, X 

s,x 


367 
384 


460 
555 


1.68 
1.84 


3.9 
5.3 


3966 
7262 


c:7o 1020 
"J ' "J 292 

1025 1827 


390 

762 3gg 


1031 i8f 
1682 iif 


278 496 
453 807 


PKS1404 


E,X 
W,X 


115 
118 


53.5 
53.7 


0.47 
0.55 


0.009 
0.011 


503 
327 


343 33I 

327 337 
316 




685 707 
DOO gg2 

487 502 

470 


4=4 468 
438 
079 384 
■^'^ 359 



Notes : 

(1) The codes for the Lobes are N — Northern, S — Southern, E — Eastern, W — Western, X — sizes from X-ray image, R — sizes from Radio image, I-EAV — Inner 
lobes in A2199, Bud — as described in Forman et al. 1 2003), Rel — Relic source (for A133 as described in .Fuiita et al. ^200Z) V 

(2) The viscosity is estimated assuming that the flow is laminar and a Reynolds number of 1000 

(3) The energy quoted here is E = PV, so the values have to be multiplied by the appropriate 7/(7 — 1). 

(4) The range on the limits on fc// from the uncertainty in the spectral index are given by the maximum values (superscript) and minimum values (subscript). 
The unceilainties from other parameters are shown in Fig.l2l 



buoyancy velocity will change with the evolution of the bubble, 
and so the timescale has uncertainties associated with it. 

These two timescales are measuring the age of the bubbles in 
two different ways, one the time for expansion to the current size 
at the sound speed of the ICM at the edge of the bubble; the other 
the time to travel to the current position at the buoyancy velocity. 
This, along with the fact that calculation for the expansion (sound 
speed) timescale uses twice the radius (see Section|2} is the most 
likely cause of the discrepancies between the timescales. From the 
geometry of the bubbles, we have assumed that the engine creating 
the bubble is not in the centre, but at one side. Therefore the part of 
the bubble which is furthest from the source, the bubbles' diameter 
away, is the part whose speed must have been less than the sound 
speed. 

In very young bubbles are the sound speed timescales less than 
the buoyancy timescales e.g. M87 and A2199 (though these seem 
to be expanding at a large distance from the central engine, and 



hence the comparatively large buoyancy timescale). The older, and 
usually larger, bubbles have buoyancy timescales which vary from 
about equal (Cygnus A) to ten times smaller than the sound speed 
timescale (A133). 

Both the buoyancy and refilling timescales depend upon the 
cluster mass within the radius the bubbles are from the centre. 
In the cases of A 133 and A4059, the masses have been esti- 
mated from a linear interpolation of the Abell radius as given in 
Reiorich & Bohrinaer ( 2002) and so the value obtained is likely to 
be accurate only to factors of two to three. In addition. ISlee et all 
(2001) also find that the travel time of the relic in A133 is longer 
than its age, a discrepancy which they propose can be resolved by 
identifying the relic not with the cD galaxy, but the one labelled 
as G in their Figure 6, which gives a travel time of around the 
age of the relic. The relic does seem to be much larger than others 
given its distance from the centre of the cluster, and so projection 
effects may also play a part. The timescales for A4059 also dif- 



© 0000 RAS, MNRAS 000, 000^00 



Particle Energies and Filling Fractions of Radio Bubbles in Cluster Cores 7 



fer by factors of up to three, and if the mass estimate is reduced 
from 1.5/1.1 X 10^"* M© for the Northern and Southern bubbles 
respectively to ~ 0.20 x lO^^M©, then the timescales become 
more comparable (x3 compared to x6 different). This is a reduc- 
tion of ~ 90%, however a linear interpolation of the mass from 
the Abell radius overestimates the masses of Perseus and A2199 
(inner bubbles) by more than a factor of ten. These new masses 
for A4059 give tbuoy = 0.71 x W jr, Uof = 2.0 x 10^ yr, 
k/fhuoy = 8.9 and fc//ref ~ 3.51 for the Northern bubble and 
tbuoy = 0.43 X 10'^ yr, t.ei = 1.7 x 10^ yr, fc//buoy = 31.9 and 
k/ frci ~ 9.92 for Southern bubble. The k/ f values do not change 
radically (~ 50% drop) as a result of this. 

The form for the synchrotron lifetime (Equation |3) is only 
valid for sufficiently large magnetic fields. Inverse Compton (IC) 
losses dominate when the energy density of the Cosmic Microwave 
Background is equal to that of the magnetic field, Ub ~ -B^/Stt, 
which corresponds to B = Bcmb . Although all the clusters anal- 
ysed are in the local Universe, it was checked whether the inferred 
magnetic field was larger than Bcmb. The lifetime determined 
from this limiting magnetic field for 1 GHz electrons (tcMs) was 
also calculated and is tabulated in the Appendix. Only A133 has a 
timescale, inferred from the sound speed, that is larger than tcMB ■ 
The magnetic field is also smaller than Bcmb , and so IC losses are 
important in this cluster. The Eastern lobe of Cygnus A has a mag- 
netic field from the refilling timescales that is less than 10% larger, 
and so IC losses may be important in this cluster as well. 

The sound speed timescale gives a lower limit on the age, and 
hence an upper limit on the magnetic field, and so the IC effects are 
likely to be important for more clusters than just the two mentioned 
above, as the actual age of the bubble is probably larger, and so the 
magnetic field smaller. 



5 DISCUSSION 

Initially it had been hoped that limits on fc// for the clusters were 
all of the same order of magnitude, however, from Fig.|2|it is ap- 
parent that they are not. The minimum possible value for fc// — 1 
is shown by the dotted line. Hydra A is the only source whose k/ f 
value falls below this line, though the uncertainties are such that 
it could have a value greater than one. A possible reason for this 
source having such a low k/ f value is that the bubbles may have 
punched their way through the ICM. Therefore the bubbles are no 
longer being contained and so are no longer in the type of equilib- 
rium that has been assumed in the calculations. 

The calculations have shown that all the lobes cannot be in 
equipartition. Although some of the limits on fc// are larger than 
the equipartition fc// values, this is due to the form of the depen- 
dence of fc// on the magnetic field (see Fig. 7 of Fabian et al. 
12002)). All the timescales derived from observational data indi- 
cate that the bubbles cannot be in equipartition. The only timescales 
that allow equipartition are the minimum and maximum from the 
dynamical constraints (Section|3}, and only for some clusters. 

The limits on fc// from the maximum and minimum ages of 
the hole as calculated from the dynamical constraints arguments in 
Section|2|are shown in Fig.|4] along with the equipartition values 
and the sound speed limits for comparison. These limits are higher 
than the ones obtained from the sound speed for many clusters. This 
implies that the assumptions used in the calculation of these limits 
on fc// may not be valid. Another indication that this may be so is 
the fact that the timescales are inconsistent with each other — the 
maximum age of the hole is frequently less than the minimum age 
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Figure 4. The value of fc// for each cluster analysed. The I, o and A sym- 
bols denote the k/ f values from the sound speed, minimum and maximum 
timescales respectively (SectionQ. The x symbol denotes the equipartition 
values. 



of the hole (see the Appendix for the values), and, on the whole, 
both are less than the sound speed values. This can be seen in Fig. 
Q where the maximum allowed age for the northern bubble in the 
Perseus cluster is less than 10^ years, whereas the values calculated 
from the sound speed and shock front arguments give 2.2 x 10^ and 
3.2 X 10^ years respectively. We therefore believe that the assump- 
tions made in the derivation and calculation of these values are not 
valid for these bubbles. The largest assumption are that the bubbles 
have been expanding sub-sonically and in a medium of constant 
pressure, which is unlikely, as there is a shock front visible in the 
Perseus X-ray images an the temperature and density of the ICM 
varies with distance from the centre of the cluster. Another assump- 
tion used is that they have not buoyantly detached from the source, 
which in some cases must have occurred (e.g. A133 and MKW3s). 

Although the fc// values calculated are upper limits, treating 
them as absolute values, it was investigated if there was any distri- 
bution of fc// with the physical parameters of the sources. The only 
parameters which came up with possible trends were the electron 
density, the 5 GHz radio power of the source and the spectral index 
of the lobes (Figs.|5||6|and0. In the plot of fc// against electron 
density there seems to be a cut-off line from the lower left to the 
upper right-hand side of the scatter plot, with all the clusters falling 
to the right of this line. There seems to be a trend with 5 GHz radio 
power, with low power sources having a high fc//, and high power 
sources having a low fc//. However the Perseus cluster is an excep- 
tion to this last trend (the points with a fc// ~ 500 at a radio power 
of ~ 3 X 10^^ W Hz"'^). Any effect from beaming has not been 
taken into account, but as the bubbles seem to be moving across the 
line of sight, it is unlikely that there are large beaming effects. A 
similar trend may be present in the plot of fc// versus a. It is also 
obvious in this plot of the large effect an uncertainty in the spectral 
index has on the uncertainty in fc//. In the above discussion the 
uncertainties on fc// values have not been taken into account, how- 
ever due to the size of the uncertainties any firm conclusions being 
drawn from the plots. We believe that there is no clear correlation 
of fc// with any physical parameter of the bubble. 

As there was no firm trend of fc// with any of the physi- 
cal parameters of the lobes or the source, other possible distribu- 
tions were looked for, and it was noticed that the fc//sound val- 
ues were clustered - the value for some clusters occurred around 
3 and others occurred around 300. The number of bubbles with 
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Figure 5. The distiibution of fc// versus electron density. The errors in 
k/ f come from the uncertainty in a. 
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Figure 8. The distribution of the numbers of bubbles with given 
logio(fc//sound) values along with the best fitting Gaussian distributions 
from least-squares analysis. 
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Figure 6. The distribution of fc// versus the radio power of the source at 5 
GHz. The errors in fc// come from the uncertainty in a. 
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Figure 7. The distribution of fc// versus a. The errors in fc// come from 
the uncertainty in a. 



a logj^Q(fc//sound) values in a given range were binned (Fig.|8}. 
The distribution appears to be bimodal, with a sample size for the 
distribution of 23. When the lower populations' values are at the 
maximum allowed by their uncertainty and the higher populations' 
values at their minimum, the two populations overlap. Therefore 
there may be a continuous distribution of fc/ /sound which we have 
not seen as our sample size is small. However, with the uncertain- 
ties that we have quoted on fc// then no conclusions can be defi- 
nite. The distribution of fc// values from all timescales was tested 
using the EMMIX^ progra m, a more recen t version of the KMM al- 
gorithm described in Ashman et aljfl994) . which uses expectation 
maximisation to perform a maximum likelihood calculation. The 
values from the sound speed argument are most likely to be de- 
scribed by a bimodal Gaussian distribution, however the buoyancy 
values are best fitted by a trimodal distribution, and the refilling 
values are all but equally likely to be fitted by two, three and four 
populations. However, all are unlikely to be described by a single 
Gaussian in log space. The binned fc// values for all of the physical 
timescale calculations are shown in Fig.|9|for comparison. 

The data from the sound speed were then fitted using least 
squares by two Gaussians, which are shown in the Fig. |8| with 
means of logio(fc//sound) = 0.44 ± 0.09 and 2.45 ± 0.06 (~ 2.3 
and ~ 280). The clustering analysis produced the following means 
log 10 (fc/ /sound) = 0.44 and = 2.51 (~ 2.8 and ~ 320), which 
are within the range from the uncertainties given by least squares 
fitting. 

In order to allow the discussion of a bimodal population of 
fc// it was investigated whether it was possible to obtain lower 
limits on fc// for the upper population. If a lower limit could be 
obtained and is larger than the upper limits of the lower population 
then the bimodality of the distribution of fc// can be discussed with 
more reliability. 

A lower limit was obtained using the ghost bubbles in Perseus, 
the only cluster where this is possible to some degree of accuracy. 
A calculation similar to that for the upper limits was used and the 
results are shown in Table |4| As there is no GHz radio emission 
from the ghost bubbles, the electrons have all aged. Therefore the 
age of the bubble must be longer than synchrotron cooling time 
of the GHz electrons. Radio emission at 330 MHz was analysed 



^ http://www.maths.uq.edu.au/~gjm/emmix/emmix.html 
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Table 4. : Lower limits on k/f from Ghost Bubbles in Perseus 



Lobe^ 


k/ /cq 


^/ /"sound 


k/f 

k/ /buoyancy 


fc/ /roflll 


''CMB 


Timescales, lO^yr 

^sound ^buoy 




West 
South 


4271 
4460 


2267 iiil 
1337 


ncn 968 
^''^ 58.6 

304 i«| 


350 gf*! 
346 ^2/ 


13.9 
13.9 


0.795 
1.59 


7.53 
7.15 


5.42 
6.27 


North 
South 


3996 
7262 


cyq 1020 

1025 iii^ 


1031 llf 
1682 


278 

453 i§I 


13.9 
13.9 


2.18 
2.15 


1.19 
1.30 


4.50 
4.90 



Notes: 

(1) The lower limits on fc// for Ghost Bubbles in the Perseus Cluster (West and South). The inner bubbles' values for the upper limits on fc// are shown for 
comparison (North and South). 

(2) The timescale for the bubble calculated from the magnetic field which produces the same energy density as that of the CMB at the redshift of the cluster 
for electrons radiating at 1 GHz 
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Figure 9. The distribution of the numbers of bubbles with given 
log]^Q(fc//) values for the sound speed, buoyancy, refilling and shock 
timescale calculations. 



to obtain the energy in electrons radiating radiating between ui = 
10 MHz and U2 = 10 GHz (see Equation^ using a steep spectral 
index of a = -2.0 ± 0.5. . 

The calculation assumes that the electron distribution is still 
described by a single spectral index over the above range of fre- 
quencies. These lower limits on fc// are above most of the upper 
limits on k/f obtained for the lower population in Figure [S] (ig- 
noring errors on the limits). There are again discrepancies in the 
timescales, mainly due to the fact that the sound speed timescale 
measures the expansion time of the bubble while the buoyancy 
timescale measures the travel time to the bubbles' current location. 
In some cases the lower limit obtained is larger than the upper limits 
for the inner lobes. However, if the ranges in the limits are included 
then there is overlap between the limits. Also, it is reasonable to 
use the limit on k/f calculated from the buoyancy timescale rather 
than the sound speed timescale as the lower limit, as these bubbles 
have detached from the centre and so are rising buoyantly through 
the ICM. Using this timescale, the lower limits on k/f are lower 
than the upper limits from the younger bubbles. 

From the observ ation that there is smooth spectral ageing 
across the inner lobes I' Fabian eta 0(2002) and no GHz emission 
from the outer ghost bubbles we are able, in the special case of the 
Perseus cluster, to obtain lower as well as upper limits on the value 
for k/f. It is constrained to be between ~ 100 and ~ 1000 with 
the assumption that there is no large change in resulting from the 



bubbles' detaching from the central engine or its age. The upper 
limits on fc// for most of the other clusters in the sample are not 
compatible with this range, which implies that there is at least a 
large spread in the values for k/f, if not a bimodal distribution. 
Although the sample size is small (23) we now go on to discuss the 
implications of a bimodal distribution. 



6 INTERPRETATION 

To obtain a high value of k/f then either fc is intrinsically high or 
/ <C 1 . If / ~ 1, then fc is high and the particles corresponding 
to the observed synchrotron radio emission do not account for all 
the particles required for there to be pressure equilibrium. Hence, 
for M84, Perseus and PKS 1404-267, all of which have high val- 
ues for k/f, the number of particles emitting synchrotron radiation 
are not sufficient for the bubble to be in pressure equilibrium. These 
bubbles, of course, may not be in pressure equilibrium, however the 
lack of particles would imply that the bubbles are collapsing, which 
assuming that they are still "powered," is unlikely. The extra parti- 
cles required could come from various sources - particles radiating 
at energies less than 10 MHz, any non-relativistic (thermal) compo- 
nent which has been swept up into the bubbles during their creation 
and also turbulence. If the jet creating the bubbles were electron- 
proton, then the protons would be undetectable from the radio emis- 
sion alone, however they would exert a pressure on the surrounding 
ICM, and so could account for the extra pressure in the bubble. This 
implies that there are two types of jets producing radio bubbles — 
heavy, electron-proton jets giving a high value for k/f, and light, 
electron-positron jets giving a low value for k/f. Both type s of jets 
have been discussed in the past — ICelotti & FabianI Jl993h present 
arguments based on the physical properties of parsec scale radio 
sourc es that radio jets are electron-proton; whereas Reynolds et aO 
il99d) argue that, at least for M87, the jets are electron-positron. 

Another effect may be due to the spectral ageing of the elec- 
trons in the radio plasma. Once the bubble plasma is no longer 
"powered" then the electrons will radiate away their energy, with 
the highest energy ones losing their energy quickest. Assuming an 
initial power law spectrum, the spectrum will steepen at the high 
frequency end. If the spectral index that has been used in the calcu- 
lations comes from the part of the spectrum that has steepened, as it 
is extrapolated back down to 10 MHz, the value for the flux emitted 
at lower energies will be an overestimate. Therefore the number of 
lower-energy particles will also be overestimated, and so pressure 
equilibrium would seem to be achieved with the "observed" par- 
ticles, so giving a low fc value, even though these particles do not 
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exist and k is high. This would imply that the low values for k/ f 
obtained for some clusters are because a spectral index that is too 
steep has been used in the calculation. MKW3s and A133 could 
be such sources. Their bubbles appear to have detached, they have 
very steep spectral indices and their k/ f values are consistent with 
one. However the radio emission from other bubbles do not have 
unusually steep spectra and hence it is unlikely that this situation 
occurs in the majority of the bubbles. 

If / is allowed to vary, then it would be possible to obtain a 
high k/f with a low k, so the particle energy factors are the same 
for all clusters, and all that varies is the volume filling factor. Us- 
ing the means of the two Gaussian distributions from Section|5| a 
volume filling fraction of ~ 0.01 for the clusters with a higher k/f 
value would mean that their k values were around 2 — 3, similar to 
the other sources. However this begs the question - why do there 
seem to be two types of bubbles - ones with / ~ 1 and ones with 
/ ~ 0.01? A bimodal distribution seems unlikely in this case, and 
therefore if the variation in fc// is due to variation on / then there 
is probably a continuum, possibly shown by fc//buoy, and the clus- 
ters that we have chosen appear to give a bimodal distribution of 
k/ f sound- If the distribution of the radio-emitting plasma were fil- 
amentary, then / may have a value ~ 0.01, and a number of radio 
lobe s do show filamentary structure, e.g. Fornax A ( Fomalont et a L, 
11989 ). M87 I Owen et al. 200C) and Cygnus A dPerlev et all 1984^ 

In Section lzTl we have calculated the effect of a simple non- 
uniform magnetic field, which approximately doubled the value of 
k/f over that obtained from a uniform field. Hence if the magnetic 
field were highly non-uniform then the k/f value may increase a 
substantial amount. As the radio emission is filamentary in some of 
the lobes which have been analysed, the magnetic field may be as 
well. 

A co mbination of the two effects may also fit the data. 
iFabian efal. (2003b,) showed that the ghost bubbles in Perseus drag 
up thermal gas (Ha) as they rise through the ICM. If this gas is 
also entrained and mixed with the relativistic plasma during the 
creation of the bubbles, then the volume filling factor of the ra- 
dio plasma would be less than one, and, there are extra particles 
present, meaning that k would be greater than one. This would also 
imply that there is some evolution of k/f with time, however plot- 
ting the k/f value against the corresponding timescale shows no 
clear correlation. However Schmidt et al. ( 2002) state that the lobes 
are essentially devoid of thermal gas, though if it were hotter than 
~ 11 keV it would not have been detected. 

The role of / could also be an explanation, i.e. whether / rep- 
resents the fraction of the volume of the bubble which is occupied 
by radio plasma, and the remainder is thermal gas, or whether / 
is the fraction of the volume filled with detectable radio emitting 
plasma, and the remainder in this case is filled with relativistic 
plasma where the electrons (and positrons) have aged, but the mag- 
netic field (and protons) have not lost any of their energy. For the 
sources with a low k/f and where the radio emission does not fully 
occupy the X-ray bubble, e.g.A2052, rather than using the flux from 
the total observes radio emission, the flux from the lower surface 
brightness region were used. Interpolating this over the whole vol- 
ume would mean that fewer particles would be inferred from the 
radio emission, and so the value for k/f would increase. 

6.1 Re-acceleration 

For there to be re-acceleration in the bubble the magnetic field has 
to be sufficiently strong. In this case the synchrotron cooling time 
of the electrons could be very short and we would not expect to 



detect GHz emission from the bubble. Nevertheless, if there is re- 
acceleration then the electrons may still be radiating, from this we 
would infer a long cooling time. Combining Equations|2|and|3]in 
the regime where the term is small enough to be ignored, then 
fc// oc tjy^^f^, and so an inferred cooling time that is long would 
lead to a limit on fc// that is lower than appropriate. 

This could explain the shape and spread of the distribution 
of the limits on k/f. If the bubbles which have a limit on k/f 
which is small have ongoing re-acceleration of electrons, then even 
though the cooling time is short, they have detectable radio emis- 
sion. Hence we have assumed that the cooling time is long, of order 
the age of the bubble, and so the calculated limit on k/f is smaller 
than if the actual cooling time had been used. In this case the bub- 
bles which have a limit on fc// which is large have comparatively 
little re-acceleration, and the limits on fc// obtained are close to 
those which have been used in the literature. 

The limits on k/f cannot be raised indefinitely - there is a 
maximum value for the limit (Section|2j. However if all the bub- 
bles whose limit on k/f is in the lower population of the bimodal 
distribution arising from the lack of strong shocks in the rims have 
their limit on k/f raised to the maximum possible, whereas the 
others remain where they are; then a multimodal distribution is no 
more likely that a single Gaussian. The increase in the magnetic 
field required to raise the limits on k/f to their maximum range 
from 1.5 (A2199, I-E) to 15 (Cygnus A, E) times that of the one 
inferred from the synchrotron cooling time with no re-acceleration. 
These fields lead to cooling times for 1 GHz electrons that are 0.54 
to 0.017 times those inferred from the age of the bubble. Hence 
many cycles of re-acceleration are required such that the electrons 
are detected, however ~ 60 cycles seems unlikely. 

There are two possibilities for how re-acceleration occurs - ei- 
ther there is a population of electrons whose number is almost fixed, 
and they are re-accelerated repeatedly many times over the age of 
the bubble so that they are still radiating today. As the radio lumi- 
nosity Lradio oc NoB^-^u^dV, then for a given luminosity, if the 
magnetic field is stronger, then there are fewer electrons present. If 
the electrons are re-accelerated, then it is likely that other particles 
present (e.g. protons & ions) will also be re-accelerated. This would 
have a minimal effect on k/f as, although there are fewer particles 
present, they have a higher energy from the re-acceleration. These 
two effects cancel each other out to some extent and so fc is about 
the same. There are no new extra particles present in the bubble in 
this case, and so / of the synchrotron emitting plasma is also about 
the same, hence there is little change in fc// and so this is unlikely 
as an explanation. 

In the other case the electrons are re-accelerated once, in a 
hotspot say, and then flow throughout the bubble where they age 
rapidly, but are continuously replaced and so there is still radio 
emission present today. Cygnus A has a radio structure which 
seems to match the latter explanation - i.e. the strongest radio emis- 
sion is at the hotspot, and then the image appears to show that the 
plasma flows away and becomes less radio bright. In this explana- 
tion there would be a population of aged electrons which reduces 
the volume filling fraction of the emitting relativistic plasma and 
increases fc as the aged electrons also exert a pressure. In this case 
we only detect ~ tsync/tbubbic of the particles present. The prob- 
lem with this model is that for the electrons to flow throughout the 
bubble before they age they have to travel at ~ 0.3 c (for a cooling 
time of ~ 10^ yr and a bubble radius of 5 kpc). 
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7 CONCLUSIONS 

From a sample of low redshift clusters with clear radio bubbles 
which are coincident with decrements in the X-ray emission, we 
have determined limits on k/ f, where k is the ratio of the total 
relativistic particle energy to that in electrons radiating between 
10 MHz to 10 GHz and / is the volume filling factor of the rel- 
ativistic plasma. We find that no bubble has simple equipartition 
between the pressures from the relativistic particles and the mag- 
netic field, k/ f was found to have no strong dependence on any 
physical parameter of the host cluster, however there seemed to be 
two populations. One set of clusters had ak/f value around 2, the 
remainder had a value of around 300. This apparent bimodality of 
the distribution of k/f could be explained in various ways. The jets 
creating the bubbles are of two types — electron-positron, which 
would have a low value for k (for / ~ 1), and electron-proton, 
where the protons are the extra particles required to maintain pres- 
sure equilibrium, but as they are unseen in the radio emission k is 
high. Spectral ageing of the radio plasma steepens the spectrum, so 
giving the impression that there are more particles present. A bi- 
modality in the volume filling fraction which could be caused by 
either a non-uniform magnetic field, or a filamentary structure in 
the lobes. If thermal plasma is entrained during the formation of 
the bubbles this would reduce the volume filling factor and provide 
extra particles, resulting in the calculated values. Variations in the 
amounts of re-acceleration may also produce the observed distribu- 
tion. 
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Table 5. Timescales and Powers 



Cluster 




,(2) 

''CMB 

lO'yr 


^sound 

lO^'yr 


^shock 

lO'^yr 


^buoy 

lO^'yr 


^refill 

lO'^yr 


lO'^yr 


^max 

lO^'yr 


sound 

1043ergs"i 


"^shock 

lO^'^ergs^i 


'^buoy 

lO^^ergs-i 


"^refill 

10"*^ergs"i 


A133 


Rel, R 


12.38 


13.9 


- 


1.99 


6.93 


1.26 


1.22 


4.74 


- 


33.1 


9.52 


A2052 


N. X 


13.21 


4.11 


- 


0.39 


1.49 


0.68 


0.26 


1.80 


- 


18.8 


4.98 




s,x 


13.21 


6.16 


- 


0.65 


2.44 


1.03 


0.43 


4.06 


- 


38.6 


10.2 


A2199 


E,R 


13.36 


4.14 


- 


1.98 


5.30 


0.50 


0.93 


3.20 


- 


6.70 


2.50 




W,R 


13.36 


5.04 


- 


1.91 


5.82 


0.58 


1.03 


4.47 


- 


11.8 


3.86 




1-E, R 


13.36 


0.30 




1.54 


1.37 


0.05 


0.24 


0.05 




0.01 


0.01 




1-W, R 


13.36 


0.45 


- 


1.26 


1.68 


0.075 


0.30 


0.12 


- 


0.04 


0.03 


A4059 


N, R 


12.65 


1.66 




0.25 


0.71 


0.15 


0.13 


0.12 




0.80 


0.29 




S,R 


12.65 


1.19 


- 


0.18 


0.74 


0.23 


0.13 


0.58 


- 


3.77 


0.93 


Centciurus 


E, R 


14.19 


0.70 




0.31 


1.02 


0.091 


0.18 


0.02 




0.05 


0.01 




E,R 


14.19 


0.81 


- 


0.46 


1.20 


0.073 


0.21 


0.01 


- 


0.03 


0.01 


Cygnus A 


E,R 


12.43 


6.53 


_ 


5.09 


11.7 


0.85 


2.06 


56.3 


_ 


72.2 


31.5 




W,R 


12.43 


6.23 




4.41 


10.1 


0.77 


1.79 


77.5 




109 


47.7 


Hydra A 


N. R 


12.57 


5.04 


- 


1.25 


3.44 


0.35 


0.61 


1.98 


- 


8.00 


2.90 




S, R 


12.57 


4.16 




2.13 


3.87 


0.28 


0.68 


1.16 




2.26 


1.24 


M84 


N, R 


14.48 


1.12 


- 


0.40 


0.74 


0.23 


0.13 


0.05 


- 


0.21 


0.08 




S, R 


14.48 


2.14 




0.56 


1.05 


0.23 


0.19 


0.03 




0.16 


0.06 


M87 


Bud, R 


14.45 


0.46 




0.55 


0.89 


0.077 


0.16 


0.20 




0.16 


0.10 


MKW3s 


S,R 


12.83 


3.11 




5.52 


5.60 


0.72 


0.99 


3.22 




1.81 


1.79 


Perseus 


N, X 


13.87 


2.18 


3.20 


1.19 


4.50 


0.36 


0.80 


5.70 


3.87 


10.4 


2.75 




s,x 


13.87 


2.15 


2.90 


1.30 


4.90 


0.36 


0.87 


7.84 


5.81 


13.0 


3.44 


PKS1404 


E,X 


13.72 


1.02 




0.40 


0.74 


0.12 


0.13 


0.03 




0.07 


0.04 




W,X 


13.72 


0.63 




0.26 


0.53 


0.14 


0.094 


0.06 




0.14 


0.07 



Notes: 

(1) The codes for the Lobes are N — No rthern, S — Southern , E — Eastern, W — Western, X — siz es from X-ray image , R — sizes from Radio image, I-EA¥ — Inner 
lobes in A2199, Bud — as described in FormanetalJ<2003h . Rel^ — Relic source as described in iFuiita et alj|2002l) . 

(2) The timescale for the bubble calculated from the magnetic field which produces the same energy density as that of the CMB at the redshift of the cluster 
for electrons radiating at 1 GHz. 

(3) The power is the PV/t work only, with — ^ not accounted for. Therefore, for a fully relativistic plasma the values for the powers need to be multiplied by 
four, and for a non-relativistic plasma, by 5/2. 

APPENDIX 

The bubble timescales and derived powers are presented in Table|5] 
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